THE EXCITATION SPECTRUM FOR WEAKLY INTERACTING 

BOSONS IN A TRAP 



PHILIP GRECH AND ROBERT SEIRINGER 

Abstract. We investigate the low-energy excitation spectrum of a Bose gas confined in a trap, 
with weak long-range repulsive interactions. In particular, we prove that the spectrum can be 
described in terms of the eigenvalues of an effective one-particle operator, as predicted by the 
Bogoliubov approximation. 



1. Introduction and main results 

1.1. Introduction. Bose-Einstein condensates of dilute atomic gases have been studied ex- 
tensively in recent years, both from an experimental and a theoretical perspective [HE]- Many 
fundamental aspects of quantum mechanics were investigated with the aid of these systems. 
One of the manifestations of their quantum behavior is superfluidity, leading to the appearance 
of quantized vortices in rotating systems [311]. This property is related to the structure of the 
low-energy excitation spectrum, via the Landau criterion [5]. Excitation spectra of atomic Bose- 
Einstein condensates have actually been measured [6] , and agreement was found with theoretical 
predictions based on the Bogoliubov approximation [7j. 

From the point of view of mathematical physics, starting with the basic underlying many-body 
Schrodinger equation, it remains a big challenge to understand many features of cold quantum 
gases I3[9]. While the validity of the Bogoliubov approximation for evaluating the ground state 
energy has been studied in several cases [10hT5] . no rigorous results on the excitation spectrum 
of many-body systems with genuine interactions among the particles are available, with the 
exception of certain exactly solvable models in one dimension [16H20] . In particular, it remains 
an open problem to verify Landau's criterion for superfluidity in interacting gases. 

In this paper, we shall prove the accuracy of the Bogoliubov approximation for the excitation 
spectrum of a trapped Bose gas, in the mean-field or Hartree limit |21|122j. where the interaction 
is weak and long-range. While the interactions among atoms in the experiments on cold gases 
are more accurately modeled as strong and short-range, effective long-range interactions can be 
achieved via application of suitable electromagnetic fields [23] • Our work generalizes the recent 
results in [23], where the validity of Bogoliubov's approximation was verified for a homogen- 
eous, translation invariant model of interacting bosons. The inhomogeneity caused by the trap 
complicates the analysis and leads to new features, due to the non-commutativity of the various 
operators appearing in the effective Bogoliubov Hamiltonian. 

1.2. Model and Main Results. We consider a system of N > 2 bosons in for general 
d > 1. The particles are confined by an external potential V ex t(x), and interact via a weak 
two-body potential, which we write as (N — l)~ 1 v(x — y). The Hamiltonian of the system reads, 
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in suitable units, 

N 1 

1=1 i<J 

with A denoting the standard Laplacian on M. d . It acts on the Hilbert space of permutation- 
symmetric square integrable functions on M> dN , as appropriate for bosons. We assume that v is a 
bounded symmetric function, which is non-negative and of positive type, i.e., has non-negative 
Fourier transform. The external potential V cx t is assumed to be locally bounded and to satisfy 
V e -xt(x) — > oo as \x\ — > oo. 

Under these assumptions on V ex t and v, the non- linear Hartree equation 

(-A + V cxt )po + (v * bo| 2 )^o = £o<A) (2) 

admits a unique strictly positive solution <po, normalized as f <Pq = 1, which is equal to the 
ground state of the corresponding Hartree energy functional. In addition, there is a complete 
set of normalized eigenfunctions {^jjieN for the Hartree operator 

H H = -A + V cxt +v*<f%. (3) 

The corresponding eigenvalues will be denoted by Eq < e± < e<i . . . . We note that ipo is 
necessarily the ground state of H^, since it is an eigenfunction that is positive. Moreover, we 
emphasize that the inequality e\ > £q is strict, since operators of the form ([3]) have a unique 
ground state [25]. This will be essential for our analysis. 
Let V denote the operator defined by the integral kernel 

V(x, y) = cp (x)v(x - y)yo{y) ■ 

As shown below, our assumptions on v imply that this defines a positive trace-class operator, 
whose trace is equal to trV = v(0) = \\ v 1 1 oo . Define also 

D := H H - e = ^(ei - e ) \ifi) (ipi\ (4) 

and let 

1 /2 

E:= (d 1 ' 2 {D + 2V)D 1 I 2 ^ . (5) 

Since V is positive and bounded, E is well-defined on the domain of D. We note that both D 
and E are, by construction, positive operators, with Dip® = Eipo = 0. The Hartree minimizer 
ifo is the only function in their kernel, all other eigenvalues of D and E are strictly positive. 

It turns out that E — D — V is a trace class operator. (We will prove this in Subsection 15.21 
below.) Let = eo < &\ < &2 < • • • denote the eigenvalues of E. Our main result concerns the 
spectrum of the Hamiltonian Hn, and reads as follows: 

Theorem 1. The ground state energy Eq(N) = inf specify equals 

E (N) = N I (IVM^I 2 + V cxt (x)ip (x) 2 )dx + ^^ I ^ {x) 2 v{x-y)<p {yfdxdy 

JR d 1 JR 2d 

- hst(D + V - E) + OiN- 1 ' 2 ) . (6) 
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Moreover, the spectrum of — Eq(N) below an energy £ is equal to finite sums of the form 

5>n 2 + O^iV-^) , (7) 

i>l 

where 6 N with Y2i>i n i — N ■ 

The error term 0(N~ 1 ^ 2 ) in ([6]) stands for an expression which is bounded by a constant times 
N~ l l 2 for large N, where the constant only depends on the interaction potential v and the gap 
£\ — Eq in the spectrum of -ffnj likewise for the error term 0{^/ 2 N~ 1 / 2 ) in (0). The dependence 
on v is relatively complicated but could in principle be computed explicitly by following our 
proof; all our bounds are quantitative. 

Our result is a manifestation of the fact that the Bogoliubov approximation becomes exact in 
the Hartree limit N — > oo. In particular, as long as £ <C N, each individual excitation energy £ 
is of the form ^j>i e « re i(l + This result is expected to be optimal in the following sense: 

if ^ <C N fails to hold then there is a non-negligible number of particles outside the condensate, 
violating a key assumption of Bogoliubov's approximation [7JEIE4] . Hence there is no reason 
why the Bogoliubov approximation should predict the correct spectrum for excitation energies 
of order N or larger. 

Theorem [1] states that the low-energy spectrum of Hn — Eq(N) is, up to small errors, equal 
to the one of the effective operator 

N 

Y,E % , E = J2e j \<p j )(<p j \, (8) 
i=i j>i 

where the subscript i in E{ stands for the action of the operator E on the i'th variable. Note 
that E is unitarily equivalent to the operator E defined in ([5]). The proof of Theorem [1] actually 
consists of constructing an explicit unitary operator that relates Hn — Eq(N) and ([8]). In other 
words, we shall bound Hjy — Eq(N) from above and below by a suitable unitary transform (cf. 
Eq. (I29p below) of ([8]), with error terms that are small in the subspace of low energy. As a 
byproduct of the proof we obtain the following corollary. 



Corollary 1. Let P 3 H be the projection onto the subspace spanned by the eigenf unctions cor- 
responding to the j lowest eigenvalues of Hn (counted with multiplicity). Similarly, let Pj^ = 

Yjk=i IV'fc) (V'fcl be the projection onto the subspace spanned by the eigenf unctions corresponding 
to the j lowest eigenvalues of 



N 



\i=l / i=l 

(k\ < /c2 < . . . ), where U is the unitary operator defined in A29\) . Then there is a constant C , 
depending only on v and e\ — e$, such that if kj + \ > kj then 

Kj + l ttj 

with || • H2 denoting the Hilbert- Schmidt norm. 
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The corollary implies, in particular, that the ground state wave function °f Hn satisfies 

^o-rf ®f =l ^f ^CN- 1 ! 2 (9) 

(for a suitable choice of the phase factor). The presence of the unitary operator IA in ^ is 
important, we do not expect that is close to ®^ 1 v , o i n an £ 2 -sense for large N. (Compare 
with Remark [U] in Section [7J) 

In addition, the corollary states that the eigenfunctions of Hn near the bottom of the spectrum 
are approximately given by applied to the eigenfunctions of ([5]), which are symmetrized 
products of the eigenfunctions ifi of Hr in ([3]). These functions can be obtained by applying a 
number, n, of raising operators a\(fi) to the N — n particle ground state, which is simply the 
product ni=i™ <Po( x i)- (Here we use the convenient Fock space notation of creation operators, 
which will be recalled in the next section.) In Subsection 15.11 we shall also calculate tfi a) \ipi)U 
(up to small error terms), and hence arrive at a convenient alternative characterization of the 
excited eigenstates of H^. (See Remark [5] in Section [71) 

Remark 1. The emergence of the effective operator E in ([5]) can also be understood as follows. 
One considers the time-dependent Hartree equation idttp = (— A + V cxt + M 2 * v)ip and looks for 
solutions of the form (p = e~ 16ot (ipo + ue~ lujt + ye lujt ) for some u > 0. Expanding to first order 
in u and y leads to the Bogoliubov-de-Gennes equations (see, e.g., [26], Eq. (5.68)) 

D + V V \(u\ fu\ (1Q) 



-v -(d + v) j v y J' V y 

The positive values which can be assumed by u are then interpreted as excitation energies. This 
is in agreement with our result: We will see below that the values for oj obtained this way are 
precisely the eigenvalues of E. (Compare with Remark H] in Section HI) 

1.3. The translation-invariant case. It is instructive to compare Theorem [T] with the trans- 
lation invariant case studied in |24j . where the Bose gas is confined to the flat unit torus T d . 
Up to an additive constant, the Hartree operator equals the Laplacian in this case, whose ei- 
genfunctions are conveniently labeled by the quantized momentum p £ (2ir r L) d , and are given 
explicitly by the plane waves <f p (x) = e ip ' x . In this basis, the operators D and V can be written 
as 

D = ^2 p 2 m (vpi 

j)G(27rZ) d 

v = "(p) fa) (fpi 

with v(p) = Jj d v(x)e~ ip ' x dx = v(—p). Since D and V commute in this case, we further have 

E= Vp 4 + 2p 2 Hp)\Vp)&p\ ■ 

Hence 



tr (D + V - E) = Y [P 2 + ~ Vp A + 2p 2 v(p) 
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and the eigenvalues of E are given by 

e p = yV + 2p 2 v(p) , 

yielding the well-known Bogoliubov spectrum of elementary excitations, which is linear in \p\ 
for small momentum. 

1.4. Short-range interactions. In Theorem [H we assumed that v(x) is a bounded function. 
If we replace v{x) by gS(x), then D + V — E will, in general, fail to be trace class (in fact, it 
is not for the above model of bosons on T d for d > 2). However, Formula ([7]) for the excitation 
spectrum still makes sense. Since all our bounds are quantitative, our proof thus shows that if 
v is allowed to depend on N in such a way that it converges to a 5-function, and v (0) increases 
with N slow enough, then the excitation spectrum is still of the form e i n i, where et are the 
non-zero eigenvalues of E in ([5]), and V is now the multiplication operator gifo(x) 2 . If v(0) 
increases too fast with N, though, our error bounds cease to be good enough to allow this 
conclusion. 

Consider now the case d = 3. If we write the interaction potential as (N — l)" 1 X^vq(Xnx) 
for some fixed, iV-independent Vq, with A^r — > oo as N — > oo, we expect that the Bogoliubov 
approximation yields the correct excitation spectrum as long as An <C N. If Ay\r ~ N, the 
scattering length of the interaction potential is of the same order as the range of the interactions. 
This corresponds to the Gross-Pitaevskii scaling [8] of a dilute gas. In this latter case, the 
scattering length becomes the physically relevant parameter quantifying the interacting strength, 
instead of J* R3 v(x)dx. Hence we expect the following to be true. 

Conjecture 1. Consider the Hamiltonian 

N 

h n p ■■= Yl (- Ai + tktfcO) + n2 ^2<N(xi - Xj )) , 

i=l i<j 

on L 2 (M. 3 )® sN , with v non-negative, bounded and integrable at infinity, and denote its ground 
state energy by E${N). The spectrum of H^ p — Eq(N) below an energy £ <C N is equal to finite 
sums of the form 

J2eini(l + o(l)) 

i>l 

for large N . Here, e, is defined as in the Hartree case with the replacements 

H H ~» H GP := -A + V cxt + 87ra ^o 
V ~> 87raoyo > 

where (fo is now the minimizer of the Gross-Pitaevskii energy functional, and ao is the zero 
energy scattering length of the interaction potential v{x). 

We expect the proof of Conjecture [1] to be more complicated than that of Theorem [TJ In 
particular, the Bogoliubov approximation would have to be modified in such a way to account 
for the detailed structure of the wave function when particles are close, which gives rise to the 
scattering length ao (instead of merely its first-order Born approximation (8ir)~ 1 J R3 v(x)dx). 
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1.5. Outline. The remainder of the paper is organized as follows. In Section [2] we estab- 
lish bounds on the number of particles outside the condensate, the iV-body Hartree operator 
YliLi D{, and their product for a low-energy state. Section [3] shows how Hn can be bounded from 
above and below by what we call the Bogoliubov Hamiltonian, which is formally close to Bogoli- 
ubov's approximate quadratic Hamiltonian on Fock space, yet is particle number conserving. 
The diagonalization of the quadratic Hamiltonian can be achieved by a Bogoliubov transforma- 
tion, which is carried out in Section [U To diagonalize the actual Bogoliubov Hamiltonian we use 
a modification thereof, which involves the estimation of various error terms (Section [5]). Finally, 
we shall complete the proof of Theorem Q] (Section [6]) and Corollary Q] (Section [7]). 

Throughout this work a multiplicative constant C in an estimate is understood to be generic: 
it can have different values on each appearance. By || • || we denote the operator or vector 
norm, depending on context; || • ||i and || • ||2 denote the trace class and Hilbert-Schmidt norms 
of operators, respectively. 

2. Bounds on the Condensate Depletion 
It is convenient to regard the iV-particle Hilbert space 

J-W := L 2 (M. d )® sN , the symmetric 
tensor product of TV" one-particle Hilbert spaces L 2 (M. d ), as a subspace of the bosonic Fock space 
T = e~ =0 J-( Ar ). The Hamiltonian H N can then be written in second quantized form as 

Hn = Yj h v a l a J + 2 fff_ ^ Vijkia}a\a k ai (11) 

i,j i,j,k,l 

where 

hij := {<pi\ - A + Vextlpj) and v ijk i := (ipi,pj\v\ip k ,<Pl) ■ 
Recall that the set {<^j}ieN denotes the orthonormal basis of eigenfunctions of Hn in ©, which 
we can all assume to be chosen real without loss of generality. The operators a\ and cij in (jlip are 
the usual creation and annihilation operators corresponding to these functions, i.e., aj := a{(fi). 

To be precise, Hn in dU) agrees with the right side of (fTT]) on the subspace F^ N \ We shall 
always work on this subspace, and use Fock space notation only for convenience. In particular, 
unless stated otherwise, all subsequent identities and inequalities involving operators on Fock 
space are understood as holding on J 7 ^) only. 

We introduce the rank-one projection P = \ipo) (y>g| and the complementary projection Q = 
1 — P. The operator that counts the number of particles outside the Hartree ground state is the 
second quantization dr(Q) of Q and will be denoted by A^ > , i.e. 

N 

i=l i 

Here and in the following, Yl denotes a sum over all nonzero indices. Another important 
quantity is the following iV-body Hartree operator, 

N 

T H := (" A * + V eAxi) + (v * (fl)(xi) - e ) = dT(D) , (12) 

i=\ 

with D defined in (gj). 

The following lemma gives simple bounds on the ground state energy of Hn, as well as on 
the expectation values of N > and Th in low-energy states. 
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Lemma 1. The ground state energy Eq(N) of Hjy satisfies the bounds 

N 1 N 

> E (N) - Nh 00 - ywoooo > 2^0000 - 2 (j\r - l) ^ ' 

Moreover, for any N -particle state ^ with (^ \ \ < S) < Nhoo + y^oooo + M> we have 

(s 1 - e ) <*| N> |tt> < <*| T H \V) < M + 2( /_ 1} t>(0) - ^ • (13) 

Recall that Eo and S\ denote the lowest two eigenvalues of the Hartree operator Hu in (J3|). 
We emphasize that e± — eo > 0. 

Proof. For the upper bound we use the trial function \N, 0, . . .) denoting a state where all 
particles occupy the ground state of the Hartree operator H^. This yields 



E (N) < {N, 0, . . .| a\ aj \N, 0, . . .) 



+ 2 (N- \\ ^ Vijkl ( N > °' ' ' a ] a l a k a i \ N , 0, • • •) 
N 

= Nh 00 + — foooo • 

For the lower bound we exploit the positive definiteness of the interaction potential v in the 
following way. With ip(x) = ^>q(x) — YliLi $( x ~ x i)> we have 

< / tp(x)v(x — y)ip(y)dxdy 
Jm?d 

N 



"^0000 



-J2(v * ipl)(xi) + — - J^vixi-Xj). (14) 



N-lj^ ™ (N-l) 



Put differently, this inequality reads 

I AT _ \ JV AT 

— — - ^2 v{xi - Xj) > — uoooo + J2{v * (fl)(xi) - — — v(0) . (15) 

i<j i=l ^ ' 

Since eo = ^oo + ^0000) we hence have 

N 

H N >Y^ {~ A i + Vext{xi) + (V * ifl)(Xi)) — V 00 00 - ^ N _ jT ^(0) 

N + 1 N 

= r H + iv/ioo + ^—^oooo - 2 ^ N _ 1 ^ (o) • 

The asserted bounds now follows, since Th > ( £ i ~~ £ o)N > ^ 0. □ 

Remark 2. The proof actually shows the operator inequality 

N + 1 N 
Tu< H N - Nh 00 — uoooo + 2 (j\f - 1) ^ 

from which (|13p readily follows. 
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In our analysis we shall also need bounds on the expectation value of the product N > Tn for 
a low-energy state. Such a bound is the content of Lemma El 

Lemma 2. Let ^> be an N -particle wave function in the spectral subspace of Hn corresponding 
to an energy E < Eq(N) + fx. Then 

(ei - eo ) I N>T H |*> < (a* - «oooo + 3i;(0)) + ^-^ v(0) - ^f^j + ±(2v(0) + yf . 

Remark 3. A slight modification of the proof yields the operator inequality 

( £l - eo)N>T u < (3u(0) - v ooo)T H + 2v(0) 2 + 2(H N - E (N)) 2 . 
Proof. We write 

A* 



(tf | N Tn |*) = (tf | |*} + (* N> (H N - E (N) 



where 5 = Y^i=i( v * ¥>o)( x i) ~~ ^ £ o - tv^T Si<j K^i — x j) + 2 + E o( N )- T ^ e second term can 
be bounded by Schwarz's inequality as 



N> (H N - E (N) -f)|*)|<f<*l (A r> ) 2 m 1/2 
For the first we use the permutation symmetry of ¥ and get 

(^|iV > 5|^) = N{V\QiS\V) . 
We split S into two parts, S = S a + Sb, where 



2<i<j 



S a := J> * <f%)(xi) - Ne - j^-j Yl ~ x ^ + f + ^o(^) 

j=2 

1 * 

S 6 := (v * <$){xx) - J^—[ Yl v ( Xl 



Using the positive definiteness of v as in (JHJ), this time for tp(x) = <j>q{x) — jrzi Si=2 ~ x i)i 
we obtain 



N 



— Y v ( x i ~ x i) > 



N 



-"^0000 



«(0) 



2<i<j i=2 

In combination with the upper bound on Eq(N) in Lemma [T] this implies that 

S a < ^ ( u (0) - ""oooo + M) • 
In particular, since 5 a commutes with Qi, we have 

N <¥| Qi5 a |*) < i 0(0) - ^oooo + y) (*l ^ > l*> • 
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To bound the contribution of Sb, we compute 

QxS b = <tf | Qi [(« * ^)(xi) - u(xi - x 2 )] |¥) 

= | [(« * ^o)0i) - «(jci - afc)] I*) 

+ (tf | QiP 2 [(« * ^o)(^l) " V(X! - X 2 )] P2 |^) 
+ I Ql^2 [(« * ^)(xi) - - X 2 )] Q 2 |*> • 

The second term on the right side of the last equation vanishes. For the first and the third, we 
use Schwarz's inequality and \(v * </3o( x i) — v{x\ — x 2 )\ < v(0) to conclude 

I (tf | QiS b I*) I < «(0) {9\ Q1Q2 \*) 1/2 + «(0) (tf | Qi I*) . 

Since 

N 2 Q1Q2 |^) < (¥| (A r> ) 2 I*) 

we have thus shown that 

I a^Th I*) < i (a* - ^0000 + 3«(o)) <*| iv> |*> 

Using Th > (ei — £o)-^ r> this implies 

<*| iV>T H I*) < (m - ^0000 + MO)) <*| iV> 1*) + ^ ( 2t, (Q)+m) 2 

4 ei — e 

The result then follows from Lemma [TJ □ 

3. The Bogoliubov Hamiltonian 

The well-known Bogoliubov approximation [7j consists of replacing the operators ao and in 
([lip by y/N, and dropping all terms higher than quadratic in the a\ and aj for i > 1. The res- 
ulting Bogoliubov Hamiltonian does not preserve particle number and is thus not suitable as an 
approximation to the full Hamiltonian fljy, as far as operator inequalities are concerned. To cir- 
cumvent this problem, we work with the following modification of the Bogoliubov Hamiltonian. 
For i > 1, we introduce the operators 

and we define the Bogoliubov Hamiltonian as 

#Bo g := E ' (£i - eo) &J&i + 5 2 % (2b\b 3 + 6(6,- + fcj&t) , (16) 

* i,3 

where Vy = uooij = i^il^l^j)- Note that this operator preserves the number of particles, hence 
we can study its restriction to J~( N \ the sector of N particles. The price to pay, as compared 
with the usual Bogoliubov Hamiltonian, is that the h,b\ do not satisfy canonical commutation 
relations, making it harder to determine the spectrum of H^og- 

In the following, we shall investigate the relation between Hn and H^og- In particular, we 
shall derive upper and lower bounds on Hn in terms of H^og, with error terms that are small 
in the low-energy sector. 
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3.1. Lower Bound. Using the positivity of the interaction potential v, a Schwarz inequality 
on p( 2 ) yields 

{P&Q + Q® P)vQ ®Q + Q® Qv{P <g> Q + Q ® P) 
> -e(P ®Q + Q® P)v(P <g>Q + Q®P)- e~ l Q ® QvQ ® Q . 

Consequently, 

v> P<S) PvP ®P + P® PvQ <g)Q + Q® QvP <g> P 

+ (1 - e)(P <g> Q + Q <g> P)v(P <8>Q + Q®P) + (1- e~ l )Q ® QvQ ® Q 

+ P ® PwP ® Q + P ® PuQ (8i P + P ® QwP ® P + Q ® PuP ® P (17) 

for any e > 0. The last term in the second line can be bounded from below by (1 — e~ 1 )v (0)Q®Q 
as long as e < 1 which we shall assume henceforth. We remark that in the case of translation 
invariance the terms in the last line vanish due to momentum conservation, but this is not the 
case here. 

In second quantized language, the bound (fT7|) implies that Hjy is bounded from below by the 
operator 



'hijctaj + y/N- 1^2' h i0 (b\ + bi) + h 00 (N - N 



(N-N>){N-N>-1) 1„ / tt 

+ ^oooo 2(N-1) + 2 ^ ij \ ij + i ' 

hi 

+ " - c ' i) N> % ~ T 0) + ^'wh 0" (jv - jv> » + {n - Ar> "'0 • (i8) 

restricted to the TV-particle sector. We note that 

i,j hi 
Since D — v(0) < —A + V cx t ~~ £ o < we can bound the last term as 

2-N> f 1 ^(^ r> ) 2 



jy - 1 ^ v u " 1 J - JV - 1 w iV 

This bound can be easily verified by investigating separately the sectors of different values of 
N > . (In particular, note that Th = on the subspace where A^ > = 0, for instance.) 
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We also have 



(^oiOi + Vij) b\bj — - 2^ ( y oiOi + Vij) 

i,j i,j 



where we have used that V as well as multiplication with v * ip^ are bounded operators with 
norm bounded by v(0). Using e = ^-oo + u oooo one verifies that 

r> , u ^ v,>^ , (N-N>)(N-N>-1) , n __ U N>(N> - l)v(0) 



e N> + h 00 (N - N>) + v 0000 ± ^7 — tt + (1 - e 



2(N - 1) v ' 2(N - 1) 

N 1 iV > (Af > - 1) 

= iV/too + yvoooo + ((!-£ >(0) + ^0000) 2 , iV _ 1 s 

The Hartree equation ([2]) implies /ijo + ^iooo = for i / 0, hence we have 

v^i £'^0 ( 6 j + 6.) + E'^= - ^) + - N> ) b i) 

i i v 

= E'^=(^ 1 -^ > ) + ( 1 -^ > ) 6 0- 

This last expression can be bounded by Schwarz's inequality: for any £ > one has 

J 1 ;^) 2 _ Ct ,( )2 NN> < fei - n>) + (i - iv>)6| 

Cv 7 ^! (iV - If/ 2 - ^ y/W^l V 71 ;i 

(1 - iV>) 2 , x9 JVAP 



C^W^T s w (iv - i) 3 /2 • 

Here we made use of 

E'koool 2 = (vol « * VoQv * Vo Ivo) < v(0) 2 . 

i 

What these computations show is that 

N 

H N > H Bog + Nh 00 + — voooo - £ £ , (19) 

where 

// 1x / x \ N>(N> - 1) 1 , x(N>) 2 

£? e = - ((1 - e^MO) + Vo0 oo) Q ;„ 7T ^ + "77 r^H + «(0)- V 



2(JV-1) N - 1 w AT 

d-iv>) 2 (ofJV > jr 2 .„ (0) i±£^iv> 

£,/F^T s y J (N - I) 3 / 2 W iV - 1 
< C (L^N" 1 + C l N~ 1 ' 2 ) (N> + 1)(T H + 1) + (iV" 1 + e + CA^ 1/2 )(T H + 1) 
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3.2. Upper bound. The upper bound on follows essentially the same lines as the lower 
bound in the previous subsection. By Schwarz's inequality 

{P®Q + Q® P)vQ ®Q + Q® Qv(P ®Q + Q®P) 

< e(P ® Q + Q ® P)v(P ®Q + Q®P) + e _1 Q ® QvQ ® Q 

and hence 

v < P ® Pt> P ® P + P <g> PvQ (8) Q + Q ® P ® P 

+ + + P)«(P ®Q + Q®P) + (1 + e -1 M0)Q ® Q 

+ P ® PwP <8> Q + P ® PuQ <g> P + P <8> QvP ®P + Q® PdP (8) P 

for any e > 0. This means that P//v is bounded from above by the expression f)18|) with e 
exchanged for — e. Using 

Yl ' fcj - EoSij) a\ aj = ^ ' ( h ij ~ £ ohj) b\bj + N N _^ Yl ' ( hi i ~ £ ° 6i ^ a ^ 



s-^> + N>T n N> 

<J2 (hij-eoSij) b\b 3 + ^Ji + v (0) " 



1,3 



N 



N -I 



we obtain 



where 



N 



Hn < H Bog + Nh 00 + — foooo + F e 



N>T H , ,3 + 2eN . ,. , , S N>(N>-1) 

+ "(0) Z , N + ((! + e ~ + ^oooo) 



(20) 



N 



N - 1 



2(N - 1) 



(1£>22 2iy> 



N 



{N- 1)3/2 



< C( (iV -1 + C" 1 ^" 1 / 2 + e" 1 ^- 1 j (7V> + 1)(T H + 1 
+ fe + CiV" 1/2 + A^ 1 ) (N> + 1) 1/2 (T H + 1) 1/2 



(21) 

Here we have again used that N > can be bounded by Th, and similarly for their square roots. 
To proceed with the analysis in Section \6.2\ it is convenient to work with the bound (I2ip on F e , 
involving only the operator (N > + 1)(Th + 1) and its square root. 

4. Symplectic DIAGONALIZATION 

In order to investigate the spectrum of the Bogoliubov Hamiltonian P/eog in (|16|) . it is useful 
to consider first the usual Bogoliubov Hamiltonian, which is the formal quadratic expression 



tfBo g = ^(V)V T 



D + V V 
V D + V 



(22) 



It is convenient to use a matrix notation where 

/ ai 
a= fl 2 

V ; 
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and T denotes transposition; e.g., a J Da^ stands for ^ -{^pi\D\ ipj)aiOj, (a)yVa stands for 

Yl'ij Vij a \ a ji e * c - The operator H-Qog is symmetric since V has real matrix elements with respect 
to the basis {<£>j}iGN- Eq. ([22]) is only a formal expression; in particular, it has an infinite ground 
state energy. It also does not preserve the particle number and hence cannot be restricted to the 
sector of N particles. Nevertheless, it serves as a useful device to motivate our analysis below 
leading to an approximate diagonalization of the actual Bogoliubov Hamiltonian H^og- 

We introduce the Segal field operators <fi = (<fii, fa, . . . ) T , 7r = (tti,tt2, . . . ) T , which are given 

C0-^(-)(iH(-')- 

They satisfy the commutation relations 

[(f>i, <f>j] = [7Tj, TTj] = , TTj] = ■ 

These remain invariant under symplectic transformations S, which satisfy 

svs = j=(_° 1 I). 

We can write 

H Bog = (<p\ir J )M 

where 

M . = l T *(D + V V \ 1JD + 2V \ 
2 V V D + V ) 2\ D ) ■ 

Here and in the following, we shall use * for the adjoint of an operator on the one-particle space 
J 7 *- 1 ) or the doubled space J 7 ^ © while we use * for the adjoint of an operator on Fock 
space. 

In order to diagonalize H^og we thus have to symplectically diagonalize M. To do so we 
introduce a real unitary operator Uq such that 

E = U£EU 

is diagonal with ordered eigenvalues, i.e. E = Yl'i e i Wi) (v 3 *! with < e\ < e2 < . . . ■ On the 
subspace QL 2 (M. d ), the operators D, E and E are invertible, and we denote their inverse by 
D' 1 , E^ 1 and E" 1 for simplicity, i.e., D^ 1 = Q(QD)~ 1 , etc. 
With 

/ D 1 / 2 \fU W E- 1 ' 2 
b ~ { D- 1 / 2 J { U J V E 1 / 2 

where A := D X I 2 E~ X I 2 and B := (A^ 1 )* , we then have S J = S* and 

S*MS = - ( E V 
2 V E J 

This corresponds to a Hamiltonian consisting of sums of independent harmonic oscillators of the 
form cf> 2 + 7if, and hence yields the desired diagonalization of H^og- 
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Remark 4. As claimed in Remark [T] in Section ll.2[ it is not difficult to see that the positive 
eigenvalues uj in (|10p are precisely the eigenvalues of E. With 



I 

Eq. (jlOp can be written as 



-i 
i 



2UTMT*il) = ujip , 

where we denote ip = (u,y) J for short. If we multiply this from the left with ITS*JT*, using 
T*I = JT* and S*J = J5 _1 , we obtain the equation 

2iITS*MST*x = u X , 

with x = T5 _1 T*^. This latter equation is simply 

E 



-E 

hence uj is indeed an eigenvalue of E, as claimed. 

The formal considerations above serve as a starting point of our analysis. Using S in ([23 
we define particle number preserving operators c = (ci , , . . . ) by 

h\ 1 / 1 i U A[/ \ / 1 1 \ / c 

6 f J 2 ^ 1 -i J V 5£/ y V -i 1 7 V c f 

1 / AJ7 + SC/o AU - BU \ ( c 

~ 2 V AU - BU AU + BUq 



4 



(24) 



Note that the operators A, B and £7o are all real, hence cj is indeed the adjoint of cj. By 
inverting 5 1 one easily obtains the inverse transformation law 

( c\_l( U^A^ + U^B-i A' 1 - U^B- 1 \(b\ 

V ct y 2 V ^o^ 4 " 1 - ^o 5-1 ^o^" 1 + us b- 1 ) ■ { 0) 

We can rewrite the Bogoliubov Hamiltonian 

^Bog = (tfy(D + v)b + h^vb + ^(tfyvtf (26) 

as a quadratic operator in these c, c T . Here, (6^) T F6 = ^ ■ Vijb\bj, etc. We insert (f24"[) into ([26]) 
and obtain 



■4 



ftog = J] 'e,4ci - ^ ' (u; (y - 1) U \ la, 

- Ij2' z v (>> c *] + [4,40 =: ( I ) + ( II ) + ( m ) > ( 2? ) 



2 



where 



y := ^i? 1 / 2 !)- 1 / 2 (D + V) D 1 /2 £ ;-l/2 + h c 
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and 



Z : 



rC/n* 



(A - B)*{D + V)(A + B) + hA + B)*V(A + B) + ^(A- B)*V(A - B) 



= -U£ [A*(D + 2V)A - B*DB-B*(D + V)A + A*(D + V)5] C7 
= i[/ * + 7)5 -B*{D + V)A] U . 

Note that Z is antisymmetric and hence 

ZijCiCj = — ZijCiCj — — Zji°i c j = ~~ ~ ^'i[ c i> c j] ■ 



1 >J 



and 



To arrive at (|27p . we have used that 

- (A + B)*(D + V)(A + B) + ^(A- B)*V(A + B) + ^(A + B)*V(A - B) 

= -[(A + B)*(D + V)(A + B) + A*VA - B*VB] 

= - [A*{D + 2V)yl + B*DB + B*(D + V)A + A*(D + V)5] 

= \ E +\ i B *( D + V) A + A*(D + V)B] 

= -E + Y 
2 



(A - B)*(D + V)(A - B) + -{A + B)*F(A - B) + -{A - B)*V(A + B) 

- [(A - B)*(D + V)(A - B) + A*VA - B*VB] 

- [A*{D + 2V)A + B*DB -B*(D + V)A -A*(D + V)B] 
-E — - [B*(D + V)A + A*(D + V)B] 

L-y. 

2 



5. Bounds on the Bogoliubov Hamiltonian 

To prove Theorem [1] we derive upper and lower bounds for the various terms (I)— (III) in (|27j) . 
This yields a bound on BB og in terms of an operator whose spectrum is explicit, as well as errors 
which are small for large N in the low-energy sector. More specifically we shall prove: 

Proposition 1. The three terms in {21^ have the following properties. There exists a unitary 
operator U : T^ N ^ — > (explicitly given in h29\) below) such that the following bounds hold 

on F^: 
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(I): For arbitrary A > we have 

> (1-X)W (y^'eiala^jU-Cil + X-^N-^^ + lXTH+l) 

^2'eicla < {1 + X)W ^'e i aJo i ^ + C(l + A- 1 )iV- 1 (iV> + l)(rH+l) 

(II): D + V — E is a trace class operator, and 
-CN- l (T H +l) 

< 2 ]T ' (u* (y - !pj U ) [ci, ct] - tr(D + F - E) + t;oooo 
< CN- l (T H +l). 

(Ill): 

-CN-'Th < Ys'Zij ([cj,ci] + [4,c]]) < CN- x Th. 



The following three subsections contain the proof of this proposition. 

5.1. Proof of Proposition [1] (I). We first refine the symplectic transformation S. By polar 
decomposition there is a unitary Wq such that A = \A*\Wq = Wo 1^4. 1- Since 

\B*\ = {A^l = (A*)" 1 (28) 



also B = \B*\W . Hence 

/ |A*|W ^o \ W 

V \B*\W U J V W 



where W = WqUq. The transformation 

W 
w 

is implementable on J 7 by a unitary W = T(VF), as it corresponds to a change of basis of the 
one-particle Hilbert space L 2 (M. d ). We define the real, bounded, and positive operator 



^logO^r 1 ) • 



Note that log |-B*| = a due to (|28p . One can show that for any t £ M the symplectic transform- 
ation 

e~ ta 
e ta 



is implemented on Fock space J- by e Xat where 
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However, it is important to note that e Xat does not preserve the particle number, and hence we 
shall instead work with e xt , where 



2 wjJ- 

1,3 

We will repeatedly need the following facts. 

Lemma 3. (i) V is a positive trace class operator. 

(ii) A — 1 and B — 1 are Hilbert- Schmidt operators. 

(iii) a is a Hilbert- Schmidt operator. 

(iv) iX : J 7 ^) — y J^( N ) is a symmetric bounded operator. 

The proof of this lemma will be given at the end of this subsection. The lemma implies that 

U := W f e x (29) 
is a particle number preserving unitary transformation on the Fock space J-, and hence we 

2 

can study its restriction to the iV-particle sector J 7 ^). With v := jj^f, G := cosh(a)W , and 
H := sinh(a)W, we define the operators di : T^ N ^ — > J 7 ^ -1 ) by 

di := a(Gipi) + ua^Hifi) . (30) 

Note that ([25]) implies 

Ci = n^—^ J2 ' {( W * cosh ( Q ))ij a A + ( W * sinn («))y a ] a o 
j 

_ a{Gipi)a\ ) a(Hipi)al 
from which we derive 



d\di = c\ci + a f (G(pi)a(G(fi) I 1 



■U _ J„. j_ „Un,„.\„fn,„A ( 1 _ a ° a o 

AT- 1 

( a o) 2fl o a o a o 



+ a(#<^)a t (#^) | 
cJq + a\Gtpi)a(Gtpi) 



(N-l) 2 N-l 
N> -2 



(aHHtpMHipi) + ||g^|| 2 ) ,\ T ' . (31) 



iV- 1 

l2 ^N>(N - N>) 
JN~T) 

The first step towards the proof of Proposition [T] (I) is the following lemma. 
Lemma 4. 

J2' e i4 d i - CN~ l T H N> < Y,' e 4 c i < Y.' e ^ di + GN ~H N> + 1 )( t h + 1 



Proof. By (j3Tj) we have 

N 



^2'eic\ci > ^2'eid\di - ^—'^2' e i a H G( Pi) a ( G< Pi) ( 32 ) 
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and similarly 



+ TW^Ty? (a\H<pi)a(H<pi) + \\Hipif) ) . (33) 



(TV -1)2 

Hence the lemma follows from the following three estimates: 

Y,'^{G<fi)a{G Vi ) <CT H , 

i 

(£%)<*(#>*) < CAr> > 

i 

^'eill^iH 2 <oo. (34) 

i 

For the proof of the first estimate note that the operator on the left side is the second 
quantization of GUqEUqG* , which we can write as 

GUqEUqG* = cosh(a)WoW * cosh(a) 

1 ( D -l/2 E l/2 D H2 E -l/2\ E ( E l/2 D -l/2 + E ~l/2 D l/2 



1 £1/2 ^-1^1/2 + £7-1/2^ E ^1/2^1 + E -l/2^ D l/2 



4 
1 

~ 4 

Since Th = dT(D) it suffices to show boundedness of the operator 

D~ 1 E 2 D~ 1 + D~ l E + ED^ 1 + 1 , 

which follows from 

\\D- 1 E 2 D- 1 \\ = ||1 + 2D~ l/2 VD~ l/2 \\ < oo . 
The second estimate in (1341) follows from the third, for which we note that 



J2'ei\\HVif = tr(HUZEU H* 



lE^WS arih(a)\\l 



t\\(D - E)D-W\\1 
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To bound the Hilbert-Schmidt norm of the operator (D — E)D 1 ' 2 , we use the integral repres- 

(i 

D - E)D- 1/2 \\ 2 = 7T 



entation x 1 / 2 = ir J °° t x l 2 I | — -^-^ ) dt, which implies that 



/■oo 

/ t 1 / 2 {(t + D 2 )- 1 -(t + E^D-^dt 
Jo 

/»oo 

= 2tt / t^it + E^D^Vit + D^dt 

JO 2 
/•oo 

<27r||y|| 2 / t 1/2 ||(t + E 2 )-^D 1/2 ||||(t + Z? 2 ) _1 |^t 
Jo 



(35) 



Using D < E and the spectral theorem one verifies that \\(t + E 2 y l D 1 / 2 \\ < ||(t + E 2 y l E l l 2 \\ < 
C(l + 1 3 / 4 )" 1 . Since also ||(t + £> 2 ) _1 || < 1/t, the integrand in the last line in (jMJ) falls off like 
i~ 5 / 4 at infinity, making the integral finite. □ 

Proposition [1] (I) is now a direct consequence of the following lemma. 

Lemma 5. For arbitrary A > we have the bounds 



W + C(l + A~ 1 )iV~ 1 (Af > + l 



Proof. We define if, := Li^ailA — di and hence, by Schwarz's inequality, 

(1 - A)W t afa i W - \- l K\Ki < d\d { < (1 + \)U^a\aU + (1 + A -1 )^//^ 

for arbitrary A > 0. We have to show that ^ f e^K^Ki < CN~ 1 (N > + 1) 2 . To simplify notation 
we define also gt := cosh(crf), ht '■= sinh(crf) and the quantity 

Kf (t) : = e- tX a(f)e tX - a(g t f) - ua\h t J) for / J_ Vo , 

which is related to Ki by Ki = Kw<pi(X)- We claim 

< CiV- 1 ^ + l) 2 (/|a 2 |/) . (36) 

Assuming this for the moment, we can use E = UqEUq as well as Wgo: 2 Wo = (log |A|) 2 to 
conclude that 

J2'eiK}Ki < CN-\N> + l) 2 J2' e i(W^i\a 2 W(pi) 

i i 

= CN- 1 ^ + l) 2 ^'( J B 1 / 2 (log|A|) 2 J E; 1 / 2 C/ ^|C/o^) 

i 

= CN^i^ + i) 2 ||£ 1/2 (iog |,4|) 2 £ 1/2 ||i . 

The claim of the lemma then follows if Hi? 1 / 2 (log | ) 2 _ZE7 1//2 1| i < oo. To see this observe that 

< - log |^| < l^l" 1 - 1 

< \A\~ 2 - 1 = E 1/2 D~ 1 E 1/2 - 1 (37) 
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which leads to 

= IK^z?" 1 - i)^ 1/2 ||l 

= \\{E-D)D~ 1 E l ' 2 \\ 2 

< \\{E - D)D- 1 / 2 \\l\\D- 1/2 E 1 l 2 \\ 2 ■ 

The claim thus follows from (|35p and boundedness of -B = D^ 1 / 2 E 1 ! 2 . 
The proof of (j36|) is a bit more elaborate. With 

[X,a(f)] = -va\af) 

we easily obtain 

K'}(t) = K a 2 f (t) -r a>af (t) 

where 

r a , v {t) := e~ tx [(1 - vJ)a{aip) - [X, u]a\lp) 
Using k/(0) = k'JQ) = 0, a second order Taylor expansion yields 



«/(*)= / ( t - s )( K ay(s)-r aJ (s))ds. 
Jo 

For any G J 7 ^ 

we introduce 

kip{t):= sup \\Kf{t) || , 
II«/II<1 

^a,V : = o SUP SU P W r a,f( s ) W) II • 
2 s<l ||a/||<l 

Note that 

sup Hk^/W l^> II = l|a|| 2 sup ||« a 2//|| a ||2(t) \ip) || < \\a\\ 2 k^(t) 

l|a/ll<l l|a/ll<l 

which yields 

K^{t) <f a ^ + \\a\\ 2 I Kip(s)ds 
Jo 

for t < 1. It follows from Gronwall's lemma (see, e.g., |27^ Thm. III. 1.1]) that 

Ml) - e" " 2 ^,^ 1 ) • 
If / G kera then Kf(t) = 0. For / £ kera 

lla/ll 

from which ([36]) follows if we can show that 

r a ^{l)<CN~ 1 l 2 \\{N> + l)\i>)\\. 
To see ([38]) we define g = af and first show that for \\g\\ < 1 

a)(ag){l - u^) 2 a(ag) < CN" 1 ^ + l) 2 
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and 



a(g)[X, v?[X, v]<J{g) < CN" 1 ^ + l) 2 



The first bound follows directly from 

t _ (2N + 3)iV> - (iV>) 2 - 5N - 1 
V ~ (N - l) 2 

and a\ag)a(ag) < \\a\\ 2 N>. To show (g0|) we write 

a(g)[X,u]1[X,u]J(g) < ( jffijffi ) <g) ' a^alaja^ a\g) 

< [ 2{ \~^ 2 +1 ) 2 a(g)\\a\\lN>(N> - l)a\g) 

< C7JV -2 [| a [| 1 1|^ || 2 iV r> (iV > - 1)(JV> + 1) 

< CN- l (N> + lf . 

To conclude the proof of f|36|) it remains to show that the inequality 

e -^ (iV > + i)V x <e ct (iV> + l) 2 
holds for t = 1. To that end we compute 

[X,N>] = ^a tJ [a taj ,N>]- va^a], N>]) 

= _J2'ai j (b i b j + b\b]). 
Taking the square of this expression yields 

[X, N>f < 2 I £ 'ctfHbj I I £'a H lW J +2 'a l3 b\b) J I £ 'a^fe, 
< Iv^v ^2 ' Oiija k iaiajalaj + 2||a|| 2 WiV > (iV > - 1) 

i,j,k,l 

(N - N>)(N - N> - 1) ( ^ , t + j. t 

<2 _ i | ^ (a )i j a\a j + 2\ 



< 2 (^xy-l) (w3w > (w > _ 1} + 4IWI2JV > + 2W3) 

+ »iMe (jy - JV>) ( y_- i y +3)+ v (Jy >-i) 

< CHallKJV* + l) 2 . 
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By Schwarz's inequality we obtain 

[X, (N> + l) 2 ] = (N> + 1)[X, N>] + [X, iV > ](iV > + 1) 

< C(N> + l) 2 

and hence it follows that 

e tx (N> + l)V tx = (N> + l) 2 + f e sX [X, (N> + l) 2 ]e~ sX ds 

Jo 

< (N> + l) 2 + C [ e sX (N> + l) 2 e- sX ds . 

Jo 

Gronwall's lemma then yields (|4ip . This completes the proof of the lemma. □ 

We conclude this section with the proof of Lemma [3j 

Proof of Lemma\^ (i): The positivity of V follows directly from the assumption that v is of 
positive type: 

(iP\V\ip)= / <p (x)(p (y)v(x - y)i()(x)i(j(y)dxdy > 0. 
Jm.™ 

In particular, the trace norm of V equals its trace, which is equal to 

trV = / (p (y) 2 v(0)dy = v(0) < oo . 

(ii): With A — 1 = (D 1 / 2 — E l / 2 )E~ 1 / 2 and the integral representation 

X V4 _ yl/4 = ^ f°° t l/4 M M ^ 

we have 

p - 1|| 2 < V2tt / t 1 / 4 !! ((t + D 2 y i - (t + E 2 y i ) E~ i / 2 \\ 2 dt 

Jo 

= 2 3 /V / t 1 ^ + &)-W 2 VD x l 2 {t + E^E^hdt 
Jo 

/•oo 

< 2 3 / 2 tt\\V\\ 2 t^Ut + D^D^WWD^it + E^E-^Wdt. (43) 
Jo 

Since D < E we can further bound 

\\D^ 2 {t + E 2 )- l E^ 2 \\ = \\E~ l / 2 (t + E 2 r l D(t + E 2 r l E^ 2 \\ l l 2 

<\\(t+E 2 r i \\<\. 

Using the spectral theorem we conclude that ||(i + D 2 )~ 1 D 1 / 2 \\ < C(l + i 3 / 4 ) -1 and hence the 
integrand in (|43p falls off like i~ 3 / 2 at infinity, making the integral finite. The estimate for 
B — 1 = D" 1 / 2 ^ 1 / 2 — L) 1 / 2 ) is obtained along the same lines. 
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(iii): We apply the integral representation logx = | (^r+i ~ x -^+t ) ^ anc ^ ^ e res °l ven t 
identity to obtain 

/■oo 

2||a|| 2 < / \\(\A*\ - (\B*\ +t)~ 1 hdt 

Jo 

poo 

< \\\A*\ - \B*\\\ 2 / 11(14*1 +t)- 1 \\\\(\B*\ +t)- 1 \\dt < oo 



o 



since \\\A*\ - \B*\\\ 2 = \\A - B\\ 2 < \\A - 1|| 2 + \\B - l\\ 2 < oo by (i). 
(iv): On J 7 ^) we have the bound 



(N 

hi h3 i,j,k,l 

<(£±?)W-U, 

which shows that X is a bounded operator. Its anti-symmetry follows directly from its definition. 

□ 

This completes the proof of part (I) of Proposition [TJ 
5.2. Proof of Proposition Q] (II). We abbreviate the symplectic transformation (|25p by 

A straightforward computation shows that 

[ci,c]} + [cj,c\] = 2 ^]yZ~~i 5i i ~ tv3T E ' ( L ii L ik ~ M jt M ik ) (a\ai + a\a k 

k,l 

We will show below that Y — E/2 and D — E are trace class, with 



r "f) = \ t ^ D + Q y ~ E ) • ( 45 ) 



Given that, we have 



tr (D + QV - E) 

l — E ' (^o (y - |) U L - M*C/ * (y - |) C/ m) ^ (a+a, + a\a k 



1,0 

N - N> 



N - 1 



Since tr QV = tr V — voooo> Proposition [T] (II) then follows if we can show that 

-CD < L*Uq [Y - ^ ) U L - M*Uq \ Y - — ) UqM < CD . (46) 



2 / " u V 2 
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We compute 

4L*C/ * (y - |) U L - 4M*C/ * (y - |) U M 

= ^D 1 ? 2 (l + D^VD' 1 ! 2 + D^Eil + D _1 / 2 yD _1 /2)j5^-i _ 2 J £>" 1 £ + h.c.) D 1 ' 2 
=: D^RD 1 ' 2 . 

Now i? is a bounded operator since DE -1 and D _1 E are bounded, which follows from 

HD.Er 1 !! 2 = \\DE~ 2 D\\ < 1 , (47) 
\\D~ l E\\ 2 = \\D~ l E 2 D' l \\ = ||1 + 2D- 1/2 VD~ 1/2 \\ < oo . (48) 



This proves (|46j) . 

We now turn to (SSI). Note that 



2{2Y - E) = B* [D + V - D 1/2 ED' 1/2 j A + A* [D + V - D~ 1/2 ED 1/2 j B . 
We claim that 

\\D 1/2 (E - D)D- 1/2 \\ 2 <oo, (49) 

\\D 1/2 {E - D)D' 1/2 + h.c.\\ l <oo, (50) 

H-D + QV- £||i < oo. (51) 

Since by Lemma [3] A — 1, — 1, are Hilbert-Schmidt and, in addition, V is trace class by 
Lemma El it follows from (j4"9]) - (f5"l~|) that 

B* (d + V - D l/2 ED~ l/2 ^j A + h.c. = D + QV - D 1/2 ED~ 1/2 + h.c. + Rest 

with ||Rest||i < oo; hence 2Y — E is trace class. Moreover, 

tr(2Y -E) = -tr ^D^ 2 (D - E)D~ 1 / 2 + h.c.) + tr 
= tr (D + — E) , 

where the first equality holds by cyclicity of the trace and the second is seen to be true by 
computing the trace in the eigenbasis of D. 
To show ([4"9"]) - ([5"T]) we compute 



D 1 / 2 {E - D)D- 1 / 2 = irD 1 / 2 / sft ((t + D 2 y l - (t + E 2 )~ l ) D~^ 2 dt 

Jo 

fOO 

= 2tt yftD(t + D 2 y l VD 1,2 {t + E 2 )- l D- l/2 dt 
Jo 

POO 

= 27t/ y/tD(t + D 2 )~ 1 V(t + D 2 )~ l dt 
Jo 

POO 

-4tt / VtD(t + D 2 )~ l VD 1,2 {t + E 2 Y l D 1,2 V{t + D 2 )~ 1 dt, (52) 
Jo 

where we applied the resolvent identity twice. The expression on the last line is trace class. This 
follows from the bound 

D(t + D^VD 1 ' 2 ^ + E^D^Vit + D 2 y l < \\D(t + D 2 )~ l f \\{t + D 2 )^ \\V\\\ , 
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where we have used that E 2 > D 2 in the second factor. The latter expression falls off like t~ 2 
for large t, making the integral finite. For the first term on the right side of (|52p . we compute 
its matrix elements. With D^ = — Eq the eigenvalues of D, 

fOO 

l dt 



'-pi 



tt / V~tD{t + D 2 y l v{t + D 2 y l i 

Jo 



VijTT 



Vt ^7 ndt = Vjn— " 

t + D 2 t + D 2 3 D l + D, 



In particular, since 



Vi4 



Di 



<\Vij\, 



the Hilbert-Schmidt property (149j) follows. Moreover 

Di 



Vi 



tJ Di + Dj 



+ j) = Vi 



i j i 



which implies (f50l) . To prove (f5T1) , one simply computes the trace of the operator in ([50]) in the 
basis of D, which leads to the conclusion that ^2i(fi\E — D\tpi) < oo. Since E — D is a positive 
operator, this implies that E — D is trace class. Since also V is trace class, this proves (151|) . 



5.3. Proof of Proposition [T] (III). Recall the notation introduced in (|44j) . A straightforward 
computation shows 



M a a k aj 



and 



k,l 



Y,'Zn ([c hCi ] + [cj,4]) = j^-j Y,' (L*ZM - M*ZL) kl (4a, + a\a k ) 



k.i 



-CD < L*ZM - M*ZL < CD . 



i.j 

Hence what we need to show is 

We observe that 

8(L*ZM - M*ZL) 

= -(B-A) (B* (D + V)A- A*(D + V) B) (A* + B*) + h.c. 
= D- 1/2 ED~ 1/2 (D + V)D 1/2 E~ 1 D 1/2 — D — v\ + h.c. 
= D 1 ' 2 ( [zr^l + D-^VD-V^DET 1 - 1 - D^VD" 1 ' 2 ] + h.c.)^ 1 / 2 . 
The operator in square brackets is bounded because of (|4"7|) and (jM]), hence the claim follows. 

6. Proof of Theorem [1] 

This section contains the proof of Theorem[TJ We split the proof into two parts, corresponding 
to the lower and upper bounds on the eigenvalues of Hjy, respectively. 
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6.1. Lower bound. By combining the estimate (|19p with Proposition [TJ we obtain the inequal- 
ity 

H N > Nh 00 + ^-^oooo + (1 - \)U ] 'wUA U + \\x{E - D — V) 

-C( (atV 1 + A^A" 1 + C 1 ^ 172 ) (N> + 1)(T H + 1) 

+ (iY- 1 +e + CA r ~ 1/2 )(T H + l)) , (53) 



which holds for any A > 0, ( > and < e < 1. Since the spectrum of Yli e i a l a i consists of 
finite sums of the form e^rii with ^ n- L < N, the desired lower bound follows directly from 
the min-max principle. In fact, for any function ^ in the spectral subspace of Hn corresponding 
to energy E < Eq(N) + £, Lemmas [U and [2] imply that 

(¥|(T H + l)|tf><C(£ + l) 

and 

(*| (N> + 1)(T H + 1) |*) < C(£ + l) 2 . 



Choosing e = 0{y£jN) = A and £ = we conclude that the spectrum of Hjy below an 
energy Eq(N) + £ is bounded from below by the corresponding spectrum of 

Nh 00 + ^-^oooo + J2' eia l ai ~ \ tl '( D + V-E)- C>(£§ iV- 1/2 ) . 

This completes the desired lower bound. 

6.2. Upper bound. A combination of (j20|) and Proposition [TJ implies that 

AT + 1 



^oooo + (1 + A)W+ f ^ 'e^ta^ W - ^tr(D + V - £) 



2 

+ C (n- 1 e~ 1 + N^X- 1 + AT" 1 + C 1 ^ 172 ) (N> + 1)(T H + 1) 
+ C (e + CA^" 172 + AT- 1 ) (AT> + 1) 1/2 (T H + 1) 1/2 , (54) 

for any A > 0, £ > and e > 0. To apply the min-max principle we need the following bound. 

Lemma 6. One has the bound 



U(N> + \){T H + l)U ] <C (y^'eiatai + lj . 



(55) 



Note that by operator monotonicity of the square root it follows immediately from Lemma [B] 
that 

U (N> + l^CZfe + l) 1/2 ^ f < C n^'e^lm + 1^ . 
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Hence we obtain from ([54 

UH N lfi < Nh 00 + ^±1« 0000 + (1 + A) 'a4 a i ~ l tT ( D + V-E) 



+ C (N- l e~ l + N~ l \~ l + N' 1 + C'N- 1 / 2 ) 'e ia \ ai + 1^ 
+ C (e + CiV~ 1/2 + ^V" 1 ) (j2 ' eia l ai + 1 j • 



(56) 



Given an eigenvalue of ^Zi&ia\ai with value £, we choose e = 0{^J^/N) = A and £ = to 
obtain iV/i 00 + ^-«oooo + £ - |tr(D + !/ — £') + 0(^ 3/2 iV" 1 / 2 ) for the right side of (p}. Hence 
the desired upper bound follows from the min-max principle. 
It remains to prove (15511 . 



Proof of Lemma If we can show that 

e x (N> + 1)(T H + l)e~ x < C{N> + 1)(T H + 1) (57) 

and 

W*DW = UqWqDWqUq < CE , (58) 

the claim follows since then 

U(N> + 1)(T H + l)W f < CW\N> + 1) 1/2 (T H + 1){N> + 1) 1/2 W 

= C(AT> + l)V2 W t( TH + i)yy(AT> + i)V2 



<C(iV> + l) ^'eio}oi + lJ 



where we have used ([57]) for the first inequality, and ([58]) for the second. 
We start with the proof of (j57|) . In fact we shall show that 

e x (N> + 1) 2 (T H + 1)V X < C{N> + 1) 2 (T H + l) 2 (59) 

from which the claim follows by operator monotonicity of the square root. We compute 

[X, {N> + 1) 2 (T H + l) 2 ] =(N> + 1)(T H + l)[X, (N> + 1)(T H + 1)] 

+ [X, (JV> + 1)(T H + 1)](JV> + 1)(T H + 1) . (60) 

With 

Ax := [X,iV>] = ^'a^fe + 6j6r 



i 



A 2 := [X,T H ] = ^ '^(ei - e ) + 
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we can bound 

[X, (N> + 1)(T H + l)] 2 = (Ai(T H + 1) + (N> + l)A 2 f 

= (A 1 (T K + 1)+A 2 (N > + 1) + [N>,A 2 ]) 2 

< C ((T H + \)A\(Tr + 1) + (N> + 1)A 2 (N> + 1) + [N>, A 2 } 2 ) . 

A\ < C\\a\\ 2 2 (N> + l) 2 

A\ < C||Da||^(iV > + l) 2 . 

[N>,A 2 ] = 2 £ ~ £ o) (b]b\ - bib, 



By ([MD we have 
and similarly 
Furthermore, since 

one checks that 



[N>,A 2 ] 2 < C\\Da\\UN> + l) 2 . 
To see that ||Z)a||2 < oo, we can proceed as in (137p and bound 

Da 2 D < D(D~ 1/2 ED- 1/2 - l) 2 D = D 1/2 (E - D)D~ 1 (E - D)D 1/2 . 
Hence we have 

\\Da\\ 2 < \\D^ 2 (E-D)D-^% 
which is finite due to (|49p . Applying Schwarz's inequality to (|60p . we have thus shown that 

[X, (N> + 1) 2 (T H + l) 2 ] < C(N> + 1) 2 (T H + l) 2 . 

We further have 

e tX( N > + \f(T n + l) 2 e~ tx = (N> + 1) 2 (T H + l) 2 

+ f e sX [X,(N> + l) 2 (T H + l) 2 ]e" sX dx 
Jo 

< {N> + 1) 2 (T H + 1) 2 

-t 



+ C f e sX (N> + l) 2 (T u + l) 2 e- sX ds 
Jo 



which by Gronwall's inequality implies (|59p . 
For the proof of (f58|) we need to show that 



WqDWq < CE 

or, equivalently, that 

D^WoE- 1 / 2 = DE-^iE-^DE-^y^E' 1 / 2 

/•oo 

= ttD t~ 1/2 (Et + D)~ l dt (61) 
Jo 
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is a bounded operator. Observe that, by (|47|) . 

\\D(Et + £>r 1 || < WDE^WWQit + DE~ 1 y 1 \\ < \\Q(t + DE^y 1 ]] . (62) 

With the aid of a Neumann expansion, one sees that the right side of (|62p can be bounded by 2t~ l 
for t > 2\\DE- 1 \\, which gives a bounded contribution to the integral in (|6ip . For t < 2\\DE \\, 
one can argue that by analyticity of the resolvent map t i— >■ (t + DE^ 1 )' 1 , as well as the fact 
that ED^ 1 is bounded, we get a uniform bound on \\Q(t + DE~ l )~ l \\. This argument does not 
yield a quantitative bound, however, since DE^ 1 is not a self-adjoint operator. To obtain an 
explicit bound, we make use of the fact that DE^ 1 — 1 is a Hilbert-Schmidt operator. In fact, 
it is even trace class, since by (j4"T|) and (f5Tj) 

\\DE~ 1 - l||i = WDE^iD - E)D~ 1 \\i < \\(D - E)D~ 1 \\i < oo . 

We shall apply the following result. 

Lemma 7 (Theorem 6.4.1 in [28]). Let A be a Hilbert-Schmidt operator. Then for z ^ o~(A) 
(the spectrum of A) 

°° \\AV k 



IK^-*)" 1 !!^ 



2 

(inW ( A)|*-*|)^Vfc! 



fc=0 

Define a to be the infimum of the spectrum of DE^ 1 on the space QJ-^ . It equals the 
infimum of the spectrum of E~ 1 ^ 2 DE~ 1 ^ 2 on that space, hence 

a=\\E 1 / 2 D- 1 E 1 / 2 \\- 1 >0. 

By Lemma [7] we thus have 

\\Q(t + DE-^W = \\Q(t + 1 + DE- 1 - I)" 1 )) 

"^ (t + a) fc +W 

V2 f\\DE- 1 
< exp 



£ + a \ t + a 

Here we have used the bound SfcLo^Vv^ — V^e^ 2 for x > (cf. p. 84 in [28]). This yields 
the desired quantitative bound, and concludes the proof of the boundedness of (|6ip . □ 

7. Consequences for Eigenvectors 
7.1. Proof of Corollary [1] We abbreviate 

oo 

H :=H N - E (N) +l=:J2hi \Xi) (Xi\ , 

i=l 

with hi < hi + \. For hj < £, it follows from (|53p and Lemmas HH2] that 

^■^^^^■(i + c^/iv) 1 / 2 ) . 

From (|56p we further deduce that < fc 3 - (l + C(£/iV) 1/2 ), and thus 

<xjl lxi> < *i (i + c(e/J>0 1/2 ) • 



.--so 
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A simple application of the min-max principle \29\ Lemma 2] then shows that if kj + \ > kj then 

jZ\{x k ^)\ 2 >J-C{i/N)^f^L. 
k,l=l j+1 j 

In other words, with P 3 K := Y%=i l^fc) (ipk\ and P h := £&=i IV>k) (ipkl 

=02 < r.(eimV* ^1=1 kl 



\P 3 K -P 3 H \\l<C{W 



kj + i kj 



This completes the proof. □ 
Remark 5. Note that the (normalized) eigenfunctions of K can be written as 

u ^(^. u \ u i\ N _ nA ...)=Y[ ^ + ^ ni tf\N-n,0,---) (63) 



where n = 5Zi>i n % 

< N, and |JV - n, 0, . . .) denotes the function ®fS x n (po £ J?( N - n ). The 
operators di are explicitly defined in (I30p . The operators are small in the low-energy subspace, 
as shown in the proof of Lemma The eigenfunctions of K (and, hence, the ones of H^) are 

thus approximately obtained by applying the raising-type operators d\ to the N — n-particle 
ground state. To explicitly estimate the difference of the functions (f63|) and 

t>l Vni ' 

however, it would be necessary to give bounds on products of powers of the operators Kj and 
d\, which are more involved than the ones used in Lemma 

Remark 6. As noted in Section 11.21 Corollary [T] implies that the ground state \I/o of is 
close, in L 2 -norm, to |iV, 0, . . .). To see the importance of the unitary operator U, one can 
calculate the matrix element 

(iV,0,...|W t |iV,0,...) = {N,0,...\e~ x \N,0,...) . (64) 

This equality follows from the fact that W leaves the Hartree ground state (po invariant. One 
readily checks that j- t (N, 0, . . .| e~ tx \N, 0, . . .) \ t=0 = 0. However, 

d 2 N 

(N, 0, . . .| e~ tx \N, 0, . . .) | t=0 = (N, 0, . . .| X 2 \N, 0, . . .) = - — _ - \\a\\ 2 , 

which is not small for large N. Hence we expect that the matrix element (|64p differs significantly 
from 1. 
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